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Abstract 

In this paper, we develop the dressing method to study the exact solutions for the vector sine- 
Gordon equation. The explicit formulas for one kink and one breather are derived. The method 
can be used to construct multi-soliton solutions. Two soliton interactions are also studied. The 
formulas for position shift of the kink and position and phase shifts of the breather are given. 
These quantities only depend on the pole positions of the dressing matrices. 


1 Introduction 


This paper is devoted to the study of an 0(n)-invariant generalisation of the sine-Gordon equation 

A (y) = “, P 2 + (<L«) = b (1) 

where the dependent variable a. = (a 1 , ■ ■ ■ , a n ) T is n-dimensional real vector held and /3 € M. Here 
and in what follows the upper index T denotes the transpose of a vector or a matrix. We use the 
notation (•, ■} for the Euclidean dot product of two vectors. 

Equation m was first appeared in [T] viewed as a reduction of the two-dimensional (){n) 
nonlinear a -model [2j. Its integrability properties were further studied afterwards. The Lax pairs 
was given in [3j and its Lagrangian formulation in [4]. Later, this equation reappeared in the study 
of connection between finite dimensional geometry, infinite dimensional geometry and integrable 
systems [5]. It was derived as the inverse how of the vector modified Korteweg-de Vries equation 


3 , __y 

'Ut — U’xxx T ~ \U > U/ U x , U — “g” > 


( 2 ) 


whose Hamiltonian, symplectic and hereditary recursion operators were naturally derived using 
the structure equation for the evolution of a curve embedded in an n-dimensional Riemannian 
manifold with constant curvature [6|. These have been recently re-derived in [7j. Besides, a partial 
classihcation of vector sine-Gordon equations using symmetry tests was done in [8]. 

Equation (jT]) is a higher-dimensional generalisation of the well-known scalar sine-Gordon equa¬ 
tion 

0 xt = sin 0. (3) 
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Indeed, it can be obtained by taking the dimension n = 1 and letting /3 = cos 0 and a 1 = sin0. 
The scalar sine-Gordon equation originates in differential geometry and has profound applications 
in physics and in life sciences (see recent review 0)- Vector generalisations of integrable equations 
have proved to be useful in applications m- They can be associated with symmetric spaces m- 
The rational dressing method was originally proposed in mm and developed in [T3j. This 
method enables one to construct multi-soliton solutions and analyse soliton interactions in detail 
using basic knowledge of Linear Algebra. In this paper, we develop the dressing method for the 
vector sine-Gordon equation dU) and show that similar to the scalar sine-Gordon equation Q there 
are two distinct types of solitons, namely kinks and breathers. One kink solution is a stationary 
wave propagating with a constant velocity. We show that a kink solution of the vector sine-Gordon 
equation can be obtained from a kink solution of d2]) by setting a = a sin (9, p = cos 9, where a is 
a constant unit length vector in W 1 . A general two kink solution of dU) cannot be obtained from 
solutions of Q, but it can be seen as a two kink solution of a vector sine-Gordon equation dU with 
n = 2. One breather solution is a localised and periodically oscillating wave moving with a constant 
velocity. One breather solution of the general O(n) invariant equation dU) can be obtained from 
a breather solution of HI) with n = 2 by an appropriate O(n) rotation. Two breathers solution 
can be obtained from the corresponding solution of dU) with n = 4, etc. Surprisingly, the effects of 
interaction, such as the displacement and a phase shift (for breathers) are exactly the same as in 
the case of the scalar sine-Gordon equation m m- Such interaction properties are naturally valid 
for the vector modified Korteweg-de Vries equation ([ 2 ]). The detailed study of soliton interactions 
for dZD when n = 2 can be found in m- 


2 Dressing method for the vector sine-Gordon equation 


In this section, we begin with the Lax representation of the vector sine-Gordon equation dU) given 
in jS], which is invariant under the reduction group Z 2 x Z 2 x Z 2 . We then study the conditions 
for the dressing matrix (assumed to be rational in spectral parameter) with the same symmetries. 
The 1-soliton solutions of dU) correspond the dressing matrix with only simple poles belonging to 
a single orbit of the reduction group. For one kink, it has two pure imaginary simple poles and for 
one breather, it has four complex simple poles. Using the dressing method, we explicitly derive one 
kink and one breather solutions starting with a trivial solution. 

The vector sine-Gordon equation © is equivalent to the compatibility condition [5] [C. A] = 0 
for two linear problems 

£T = 0, AT = 0, (4) 

where 


C = D x — XJ — U and A = D t + \~ L V, 


( 5 ) 


and 


0 1 0 
J = I -1 0 0 


, U = 


0 0 (L 


0 0 0 T 

0 0 -c%/P 

0 a x /P 0 n 


0 (3 a 1 

,v=\ -p 0 o 1 

—a 0 0„ 


( 6 ) 


where 0 is n-dimensional zero column vector and 0 n is the n x n zero matrix. Without causing 
confusion, we sometimes simply write 0 instead. 
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The Lax operators £ and A are invariant under the (reduction) group of automorphisms gen¬ 
erated by the following three transformations: the first one is 

1 : £(A) —»■ —A (A), (7) 

where A(A) is the adjoint operator defined by A (A) = — D x — A J T — U T . The invariance under 
this transformation implies the matrices J and U are skew-symmetric. The second one is 

r :£(A)^£(A), (8) 

where £(A) means its complex conjugate. The invariance under this transformation reflects that 
the entries of matrices U and V are real. The last one is called Cartan involution 

s : £(A) -> Q£(-A)Q, (9) 

where Q = diag(—1,1,... , 1), which leads to the reduction to the symmetric space. 

These three commuting transformations satisfy 

2 2 2-i 

l = r = s = id 

and therefore generate the group Z 2 x Z 2 x Z 2 . Indeed, the operator A is also invariant under it, 
that is, 

i{A{X))=A{X), r(A(X)) = A(X), s(A{\)) = A{X). (10) 

Thus we say the Lax representation of JT|) is invariant under the reduction group |13) 16j ]T7] 
Z 2 x 7 L 2 x 21 2 ■ 

In what follows, we use the method of rational dressing mmm to construct new exact 
solutions of dU) starting from an exact solution c?o, /3o■ Let us denote by Uq,Vq the matrices U. V 
in which a, /3 are replaced by the exact solution do, Po °f ©■ The corresponding overdetermined 
linear system 


A)'Ll) — (D x — XJ — C/o)'Lo — 0 Alo'I'o — (D t + A 1 Vd) v I r o — 0 (11) 

has a fundamental solution To (A ,x,t) invariant under transformations ([7]) (El) • Following [2l 112] 
we shall assume that the fundamental solution T(A ,x,t) f° r the new (“dressed”) linear problems 

£T = (D x -XJ- t/)T = 0 AV = (D t + A -1 V)T = 0 (12) 

is of the form 

T = T(A)T 0 , detT/0, (13) 

where the dressing matrix T(A) is assumed to be rational in the spectral parameter A and to be 
invariant with respect to symmetries 

T(A)- 1 = T t (A), T(A) = T(A), QT(-A)Q = T(A). (14) 

Conditions (1141) guarantee that the corresponding Lax operators £ and A are invariant with respect 
to transformations (ED-®. 
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It follows from (fTITl . (fl2l) and (fT3|) that 


$(D X - A J - Uo)®- 1 = -A J - U; (15) 

$(A + = A'V. (16) 

These equations enable us to specify the form of the dressing matrix and construct the corre¬ 
sponding “dressed” solution a, /3 of the vector sine-Gordon equation (HD- 

Let us consider the most trivial case when the dressing matrix $ does not depend on the spectral 
parameter A. In this case the dressing results in a point transformation (0(n) rotation) of the initial 
solution do- 

Proposition 1. Assume that & is a A independent dressing matrix for the vector sine-Gordon 
equation W- If it is invariant with respect to symmetries ( 0 , then 

/ 1 0 0 T \ 

T = ± [ 0 1 0 T , (17) 

\ o o n ) 

where 11 G 0(n,M) is a constant (x. t- independent) matrix. The corresponding “dressed” solution 
is a = Q do, (3 = Po■ 

Proof. Under the assumption that the dressing matrix <f> is independent of the spectral parameter 
A, it follows from (fT5j) and (fTUl) that 

[J, $] = 0 ; <& x = IM> - W 0 ; = 0; V<f> - 4>U 0 = 0. (18) 

It follows from the first condition [J, 4?] = 0 that the matrix <f> is of the form 

/ mu mi2 0 T \ 

4> = —mi2 mn 0 T , 

\ o on) 

where 11 is an n x n matrix. Substituting it into the last condition in (1181) . we get mi 2 = 0 and 

(/3 — Po)mi\ = 0; a T Sl = mnaJ; mua = Hao- 

This implies that (3 = f3 o since det 4> ^ 0 and nn T = rri \j I n . where I n is the n x n identity matrix. 
The rest two conditions in (1181) imply that matrix <3? is independent of x and t. Finally, since 4> 
satisfies (ED, we have that matrix n is real and m(, = 1 , i.e. mu = ±1. □ 

A A-dependent dressing matrix 4>(A), which is invariant with respect to symmetries (1141) has 
poles at the orbits of the reduction group. Simplest “one soliton” dressing correspond to the cases 
when matrix 4>(A) has only simple poles belonging to a single orbit. There are two non-trivial cases: 

(i) Matrix < 1>(A) has two pure imaginary poles at hiv, v ^ 0, i/gl. 

(ii) Matrix 4>(A) has four complex poles at ±/j, ±/i, /jl = 7 + iu, qcu 7 ^ 0, 7 , cj G M. 
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These cases correspond to the “kink” and “breather” solutions respectively. We will show that the 
case of two real poles at ±/r € M leads to a trivial dressing. 

A more general case in which the poles belong to a finite union of orbits corresponds to a multi- 
soliton solution. Using this approach we can construct explicitly multi-kink-breather solutions and 
analyse the result of kink-kink, breather-breather and kink-breather collisions. 

The explicit forms of matrix 4>(A) corresponding the above two cases and satisfying the last two 
conditions in (1141) are 


(i) 


n 


$(A) = C + —— 
A — fi 

4>(A ) = C + — 

A — n 


QAQ 
A + /./ 
QAQ 
A + /r 


A — —QAQ, 

A QAQ 
A — /r A + /u 


H = iv; 


fj, = 7 + iui, 


where C = QCQ is a real matrix of size (n + 2) x (n + 2) independent of A. 

We now derive the conditions on matrices A and C such that 4>(A) is an invariant dressing 
matrix, i.e. satisfying conditions (|14l) . First we have 

Proposition 2 . If 4>(A) is a dressing matrix satisfying 4>(A)4? r (A) = I n + 2 , then matrix C is a 
constant matrix of the same form as 

Proof. All identities must be satisfied for all values of A. For (|15l) . the linear terms in A as A —>• oo 
leads to [ J, C] = 0. So matrix C is of the form 

/ cn Cl 2 0 r \ 

C = j — C12 Cn 0 r I , 

\ o o n J 


where Q is an n x n matrix. Substituting it into the last condition C = QCQ, we get C 12 = 0. 

The constant terms of (fl5l) in A as A —» oo leads to 

C X = UC-CU 0 -[B,J], (19) 

where B is the coefficient at A _1 in the expansion <1>(A) = C + B\ -1 + • • • as A — > oo. Substituting 
the form of C into it, we get C x = 0 implying that matrix C is independent of x. For (USD and 
4>(A)<I )i (A) = I n + 2 , taking the limit A — >• oo, we get Ct = 0 and CC 1 = I n +2■ Thus we proved the 
statement. □ 


Note that a composition of two dressing matrices is also a dressing matrix. Therefore, without 
loss of generality, from now on, we take C = I n + 2 in the dressing matrix 41(A). 

Proposition 3. The dressing matrix satisfies 4>(A)4> T (A) = I n+ 2 if and only if 

AA t = 0; PA t + AP t = 0, 

where 


( 20 ) 

( 21 ) 


5 









Proof. The product <h(A)^> T (A) is a rational matrix function of A and at A = oo it is equal to 
the unit matrix. Now it is sufficient to show that the rest of conditions in (1201) are equivalent to 
vanishing all other poles of the product. The first condition of (1201) is equivalent to vanishing of 
the second order pole at A = //. All other second order poles of the product vanishes due to the 
reduction group symmetry. The residue at A = p vanishes if and only if A satisfies the second 
condition of (l 20 l) . □ 

Therefore, the matrix A is degenerate. Here we shall study solution of rank one, that is, 
rank(A) = 1 and thus represent the matrix A by a bi-vector 

A = b)(a, (a = (p, q, a 1 ,..., a n ), b) = (b 1 ,..., b n+2 ) T , (22) 

where p,q,a l ,b l € C. We use different notations for the first two component of the vector (a in 
order to emphasise that they play a particular role in the solutions. 

Under the assumption that A is a bi-vector, we rewrite conditions (1201) for A into conditions 
for vectors a) and b). 

Proposition 4. Let A = 6 )(a / 0 be a bi-vector. The dressing matrix satisfies < h(A) < l ):r (A) = /„+ 2 
if and only if 


(aa) = 0; Pa) = 0, (23) 

where P is given by H21\) . 

Proof. The first equation of (1231) immediately follows from the first equation of (1201) . For matrix 
A = b) (a the second condition of (1201) is equivalent to 

pa t = 0 . 

Indeed, the first term of the sum in the second equation of (1201) is a bi-vector with image and 
co-image spaces spanned by b ) and (b respectively and thus it is equal to 5b) (b , that is, 

PA t = 5b) (b 

for some 5 € C. The second term of the sum in the second equation of (1201) is just the transposition 
of the first one and thus we get 

PA t + AP t = 25b)(b = 0. 

Thus <5 = 0, that is, PA T = 0 implying Pa) = 0. □ 

We now investigate the conditions (1151) and (fTUl) that <F(A) must satisfy in order to be a dressing 
matrix. Notice that they are rational matrix functions of A. We compare the residues of all poles 
of both sides. For (1151) . we know that the linear terms in A vanishes as A —>• 00 from Proposition [ 2 ] 
and its constant terms in A as A —» 00 GSD becomes 

U = U 0 + [B,J\, (24) 

where B is defined by <f>(A) = I n+ 2 + B A ” 1 + • • • as A —>• 00 . We now compute residue of (fT51) at 
A = p. This leads to 

- AJA t + A(D X -pj - U 0 )P T + P(D X - pJ - Uq)A t = 0, (25) 
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where A = b) (a ^ 0 satisfying ([23]). In (1251) the first term vanishes since A is a bi-vector and J is 
a skew-symmetric matrix. It follows from vanishing of the second order pole at A = /i that 

A(D X -fiJ- U q )A t = 0. (26) 

Thus we can see that 

£ 0 a) = (D x - nJ - U 0 )a) = 0 (27) 

is a solution of (1251) and (1261) due to (|23l) . In fact, the general solution should be Cqcl) = ya), where 
7 is an arbitrary scalar function of x and t. Without any loss of generality we set 7 = 0. Indeed, 
the 7 can be removed by the rescaling a) to exp (f ydx)a) and it will be shown that the matrix A 
is not affected by the scaling of a). 

We carry out the similar analysis for (11611 and obtain that 

V = d>(0)Vo < I) T (0); (residue at A = 0) (28) 

A(Dt + fi~ 1 Vo)A T = 0; (double pole at A = //) (29) 

— h~ 2 AVqA t + A(Dt + n~ 1 Vo)P T + P(D t + h~ 1 Vq)A t = 0. (residue at A = n) (30) 

For A = b) (a, it is obviously AVqA t = 0 since Vo is a skew-symmetric matrix. In the same way as 
we get (1271) we can see that 

(A + lT l Vo)a) = 0 (31) 

is a solution of (17CT1) and (12U1) . Thus for given constant vector a\) satisfying (oyai) = 0, we obtain 
that 

a) = T 0 (x,t,/i)ai). (32) 

To construct the exact solutions for the vector sine-Gordon equation m using (1281) or ([Mil , it is 
required to determine the vector b) in the bi-vector A using (12311 . The latter depends on the choice 
of the forms for ^(A). We will determine the vector b) for the kink and breather solutions in the 
following sections. 

In the following two sections 12.11 and 12.21 we construct the exact solutions starting with the 
trivial solution /3q = 1 and ao = 0 for the equation JT]). Then Uq = 0 and Vo = J■ It is easy to see 
that in this case the fundamental solution for m is 

( cos 9 sin 9 0 \ 

— sin0 cos 9 0 J , 9 = Xx — X~ 1 t. (33) 

0 0 I ) 

Matrix To(x,t, A) obviously satisfies the reduction group symmetry conditions (1171) . 


2.1 A kink solution 


Let us assume that the dressing matrix is regular at A = 00 and has only two simple poles at points 
A = ±(i ^ 0. The set of poles is an orbit of the reduction group if /r is either real or pure imaginary. 
Then, without any loss of generality, it can be written in the form 


$(A) = /„+2 + 


A 

A — n 


QAQ 
A + n 


(34) 


For the bi-vector A = b)(a ^ 0, we now determine the vector b) using Proposition [4j 
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Proposition 5. Let A = b) (a ^ 0 be a bi-vector . The dressing matrix Ij34\ ) satisfies < F(A)<I> T (A) = 
I n +2 if and only if 


(aa) = 0 ; 


j \ _ 2/jQa) 

(aQa) 


(35) 


Proof. The first equation of (1351) has been proved in Proposition [4j For matrix A = b) (a the second 
condition of (1231) is 


Qb)(aQ . Qb)(aQa) 

( 4+2 -x- )a) = 0 , that is, o) =---. 

Zp Zp 

This leads to the second equation of (f35j) since Q 2 = 4 + 2 • □ 

Remark 1. It follows from this proposition that the bi-vector A = b)(a is parametrised by a subspace 
spanned by vector a) rather than vector a) itself. Indeed, scaling a) e+ fia) does not change matrix 

A. 


Remark 2. As we mention above, there are two possibilities for a two points orbit, namely p = 
and p = ±iu, v € M. In the first case to ensure that the solutions are real we require that the 
matrix A is real and consequently the vector a) is proportional to a real vector a) = fid), where 
a) € M n+2 . For real a) it follows from the condition (aa) = 0 that a) = 0 and so A = 0. Thus the 
case p = ±+, v € R leads to a trivial result. 

We shall see that the case p = ±in, v e R and A = —QAQ yields a nontrivial solution. 

To get the kink solution similar to the one for the classical sine-Gordon equation, we take p to 
be pure imaginary, that is, p = iu, where v € M as we stated in Remark O 

Proposition 6. A kink solution of the vector sine-Gordon system m on a trivial background 
( 7 ?o = 1 , Oq = given by 

13 = 1 - 2(p\ - q\)(pi coshp + qi sinhp) -2 , (36) 

a = 2ai(pi sinhp + qi cosh p)(p\ cosh p + q\ sinhp) -2 , i = l,...,n. (37) 

Here p = vx + o~ 1 t with v € M and (aq = (ipi, qi,a[), where pi,qi € M and ai e 1" is a constant 
vector satisfying (aia\) = 0. 

Proof. When p = iu, the fundamental solution Tq given by (1331) becomes 


'&o(x,t,iu) 


cosh p 
—i sinh p 
0 


i sinh p 
cosh p 
0 



p = ux + v 1 t.. 


(38) 


It follows from (1321) that 

a) = *„(*, iv)a\) = (i(p\ cosh p + q\ sinh p),p\ sinhp + q\ cosh p, af ) T . 


Thus 


(aQa) = 2(pi cosh p + q\ sinh p) 2 . 





It follows from Proposition [5] that 


b) = 


iv 


T\T 


(pi cosh p + qi sinh p) 


cosh p + qi sinh p),pi sinhp + q± coshp, a() 


We can now write down the matrix A. Further we use formula (1281) with Vo = J and get ()37D and 
(1551) for a and /3 as in the statement. □ 


One kink solution (1371) belongs to a one dimensional subspace (it is proportional to the vector ai) 
and thus it can be easily reconstructed from the kink solution of the scalar sine-Gordon equation. 
Indeed, if 6 is a kink solution of equation 0 x t = sin@, then 

(3 = cos 6 , a = —sin 6 
l a i| 

is the corresponding solution of the vector sine-Gordon equation. 

Later we shall show that a general two kink solution belongs to two dimensional subspace 
Span K (ai,a 2 ) and it cannot be obtained from solutions of the scalar sine-Gordon equation. In 
particular, for two-kink solutions we compute the position shifts after collision. 

For one kink solution in Proposition [ 6 j we rewrite their denominator in the following form: 

(pi cosh p + qi sinhp ) 2 = ( — gl \ (cosh (2v(x + v~ 2 t - x 0 )) + 1), = 77 — In I —-—I . (39) 

V 2 ) 2z/ pi + qi 

Thus the kink is moving on the (x, t)-plane along the straight line and its position at time t is 

x = — v~ 2 t + Xq, 

which is only dependent on the pole p = iv and the first two components of vector ai). The width 
of the kink is of the order v~ l . 


2.2 A breather solution 


A breather solution corresponds to the only simple poles at points of a generic orbit of the reduction 
group. Assuming that the dressing matrix has a simple pole at a point p = 7 + iw , p / ±p, where 
we use p to denote the complex conjugate of p, we have 


$(A) = I n +2 + 


A 

A — p 


QAQ | A 
A T p A — p 


QAQ 
A + p 


(40) 


Assume that A = b) (a ^ 0. We determine the vector b) using Proposition [3] and construct the 
solution on a trivial background (/3 q = 1 , ckq = 0 ). 


Proposition 7. Let A = b) (a ^ 0 be a bi-vector, where (a and (b are two complex vectors given 
by (a = (p,q,a l , ■ • • , a n ) and (b = (b 1 , b 2 , ■ ■ ■ , b n+2 ) respectively. The dressing matrix (JW satisfies 


^(A )^ 2 (A) = I n - 1-2 if and only if 

(aa) = 0, (41) 

(42) 

b 2 = _ 2 ^ + 2 uj 2 \p\ 2 ) - A'yuplnfipq)) , (43) 

b k +2 _ ^a k (jl 2 (aa) + 2u 2 \p\ 2 ) — 4'yuplm(pa k )^ , k = l,...,n, (44) 
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where 


A = \p\ 2 (aa) — 2u 2 \p\ 2 


(45) 


Proof. Equation (14111 has been proved in Proposition |4j For matrix A = b)(a the second condition 
of (1231) is 

„ x _, T QAQ , A QAQ , x _ n 
Pa) — (I n +2 — b 7 r — , _ )a) — 0 , 


that is, 


2p n — n p + p- 
Qb)(aQa) b)(aa) Qb)(aQa) 


a ^ U/ \u-u-7 _ ^ 

2 p p — p p + p 

We solve this linear system of equations (14611 for unknowns b l ,, b n+2 by taking into account (14111 
and obtain the formulas in the statement. □ 


Proposition 8. A breather solution of the vector sine-Gordon system m on a trivial background 
(Po = 1, cfo = 0 ) is given by 

P = 1 - -Jr- ((««>M 2 -2(Im(pg)) 2 ) , (47) 

°k _ _)_ 2uj 2 \p\ 2 )qa k ) — 4yw Im (qp) Im(pa fc )^ , k = 1 ,..., n , (48) 

where A is given by m and the components of the vectors (a = (p, q, a 1 ,..., a n ) are the following 
functions of x,t given by 

p = cos(cj) + irj)pi + sm(<f> + ir))qi, 

q = cos(0 + irf)q\ — sin(<(> + irj)pi, (49) 

a k = a\, k = 1 ,..., n, 

where <f> = yx — V = u;.x+ vector (oi = (pi, gi, a \,..., af) is an arbitrary complex 

constant vector satisfying the condition (oyai) = 0. 

Proof. When p = 7 + iu, the fundamental solution 'ho given by (1441) becomes 

( cos((j> + irj) * sin(</> + ip) 0 \ 

—i sin ((j) + ip) cos (4> + ip) 0 , (50) 

0 0 I ) 


where cj> = 7 (x — (y 2 + a; 2 ) 1 t) 
It follows from (12811 that 


77 = co {x + (y 2 + u 2 ) 1 t ). Using (l32ll we get (H9l) for a). 
V = $(0)J$ T (0), 


where 4*(A) is defined in (14011 . Making substitution of A = b) (a in (14011 and taking into account ([6]) 
we obtain 


P = Vi, 2 = |“|4(H 2 - 4 Re(//p6 1 ))( \p\ 2 - 4R e{pqb 2 ))- (51) 

oik = Vi,fc +2 = -t-tiGmI 2 - 4Re(/lp6 1 ))Re(/lg6 fc+2 ), /c = l,---,n. (52) 
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From (14211 in Proposition [7] it follows that 


R e(Hrb') = Re Up = 0. 

Thus ( \fi \ 2 — 4R e(flpb 1 )) = \p\ 2 , and therefore 

P = 1 - An Re (Mb 2 ), 

Oik = --An Re{fiqb k+2 ), 

m 

The substitution of (H3l) and (Hill in the above equations leads to (H71) and (HHl) . □ 

The breather solution (H71) . (HHl) represents a periodically oscillating localized wave of a charac¬ 
teristic width ~ w -1 moving with the constant velocity (j 2 +o; 2 ) _1 . Thus the width of the wave, 
its speed and frequency of oscillations depend on the position of complex pole p only. The location 
of the wave and the phase of its oscillations are determined by pi,qi, i.e. the first two entries of 
the vector ai). These can be found using the denominator A defined by (1451) as follows: 

A = |/i| 2 ((aiai) - |pi| 2 - |gi| 2 ) + 7 2 (|pi| 2 + |<?i| 2 ) cosh(2 rj) - 2y 2 Im(pi<?i) sinh(2 rj) + 

+uj 2 (\qi\ 2 - |pi| 2 ) cos(20) — 2w 2 Re(pigi) sin(2(/)) 
n 

= (7 2 + c j 2 ) ^2 I°11 2 + \pl + <?i I (t 2 cosh(2 rj - 2 r) 0 ) - u 2 cos(2 4> - <£ 0 )) , 

k= 1 


k = 1 ,..., n . 


(53) 

(54) 


where 



Pi ~ qii 
Pi + qii 


and 4 >o satisfies cos 4 >q 


bil 2 -kil 2 

[pf+iil 


and sin 4>o 


2 Rc(p i gi) . ypl 1U s we have 
\pI+iI\ 


(55) 


4 > o = arg 


( pi + 

\pi - qii) 


(56) 


Thus the breather is moving on the (x, f)-plane along the straight line and its position at time t is 
given by 


x = -(7 2 + w 2 ) 4 + xo, x 0 = uj 1 i] o, 

which is only dependent on the position of the pole // and the first two components of vector ai). 
We define the phase of this breather as 4> o- In Figure [TJ we plot one breather solution when n = 2 
for ai and <37 with fi = 2 + and (ai = (i, 2,3 i, \/6) at time t = 3. The position of this breather 
is at x ~ —2.936 and its phase = 0. 

The breather solution (0ZD, (USD also depends on the constant complex n dimensional vector 
ai = (aj,... , a'l) T € C n , which is the last n components of the vector ai) = (p±,qi,a \,... , a™ ) T . 
In general, the vector a of the solution (1481) belongs to the two dimensional subspace VF(ai) = 
Span K (Reai,Imai) C M n . If the subspace W(sl\) is one dimensional (Reai ~ Imai) then the 
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obtained one breather solution represents the well known breather solution of the scalar sine- 
Gordon equation (with 6 x t = sin (9 : /3 = cos 9, a = sin 6 ), which is a subsystem of the vector 
sine-Gordon equation. Solution (H8l) corresponding to a two-dimensional subspace W (a) cannot be 
obtained from a solution of the scalar sine-Gordon equation. However, it can be obtained from the 
breather solution of the two-dimensional vector sine-Gordon equation. 

Similar to the kink case (see Remark [1]), the breather solution is parametrised by a point on 
the complex Grassmanian Gi'i^+n = CP n+1 rather than the complex vector a) itself. 


3 Multi-soliton solutions and soliton interactions 


The vector sine-Gordon equation has two types of one soliton solutions, namely the kink and 
breather solutions discussed in sections [2TT1 and m The set of poles of the dressing matrix <h(A) 
for multi-soliton solutions is a union of a finite number of the kink and breather type orbits. We can 
construct the dressing matrix recursively by presenting it as a composition of elementary dressing 
factors of the kink and breather type. Another approach is based on a representation of the dressing 
matrix in the form of partial fractions with respect to the spectral parameter A. The latter approach 
is more conventional, but leads to a big system of linear algebraic equations, the size of which is 
determined by the number of solitons in the multi-soliton solution. 

In order to construct a general multi-soliton dressing of the vector sine-Gordon equation JT]) 
with m\ kinks and m 2 breathers we need to choose the following data: 

• a set of distinct positive real numbers {z 7 ,}^, 

• a set of m\ vectors a s ), s = 1,..., mi of the form a s ) = ( ip s , q s ,al,..., a ™) T , such that 

(a s a s ) = —p 2 + q 2 s + (a «) 2 + • • • + (a ”) 2 = 0 , p s , q s ,a],... ,a™ € M 
and the real “sub-vector” (a|,,..., a”) is non-zero. 

• a set of m 2 distinct complex numbers {p s = 7 s + iu s \ ,j s > 0 , oj s > 

• a set of m 2 complex vectors a s ) = (p s ,q s ,al,... ,a ™) T , s = mi + l,...,mi + m 2 and 

p s ,q s ,al,... ,a™ € C such that ( a s a s ) = 0 and “sub-vectors” ( p s ,q s ) and are 

both non-zero. 


Let us first construct the ?ni-kink solution. We start from a trivial solution of the vector sine- 
Gordon equation = 1, ao = 0, so that 

Uq = 0, Vq = J, Tq(A, x, t) = exp(AxJ — A _1 t J). 


We denote the kink dressing matrix constructed in the Section 12.11 as 

1 Qa)(a 1 a)(aQ 


<3?(A, p, a)) — I n +2 + 


(A - p) (ciQa ) (A + p) (aQa) 


Now the multi-kink solution of the vector sine-Gordon equation and the corresponding the funda¬ 
mental solution for the associated linear problems can be found recursively 


V s = $(0, iv s , ^ s -i(iv s , x,t)a s ))V s -i® T {0,iv s , ^ s -i(iv s , x,t)a s }); 
’MA ,x,t) = ^(A ,iv s , T s _i(w s ,x,t)a s ))T s _i(A,x,f); 


s = 1,2,...,mi. 
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Having constructed the m-i-kink solution V mi and the corresponding fundamental solution of the 
linear problem (A, x, t ) we can add m2 breathers using the breather dressing matrix constructed 

in Section 12.21 (see Proposition [7]) denoted as 




I-n+2 + 


A 

A — n 


QAQ A 
A T /i A - fi 


QAQ 
A + fl 


in a similar recursive way 


V s = <f >(0 ,ii s ,^s-i{^ s ,x,t)a s ))V s -i^ T ( 0 1 ii s ,^ s - 1 (/i s ,x,t)a s ))-, 

^s(A ,x,t) = ®(\,[i s ,'S> s -i(tJ-s,x,t)as))'i> s -i{\,x,t)-, 


s = mi + 1,mi + m2. 


Moreover, we could change the order of dressings by creating m2-breather solutions and then adding 
mi kinks, or even make dressings in an arbitrary order. Such recursive approach is useful for the 
study of the effects of soliton interactions as shown in the following sections. 


3.1 Interaction of two kinks 

In this section, we study the interaction of two kinks. A two-kink solution corresponding to the 
sets 

-j>i,ai) = (ip^q^al,... ,a r Q T }, {v2,a 2 ) = (ip2, 92 , < 4 , ..., a£) T }, 0 < v x < v 2 

represents two kinks moving to the left with speeds zq“ 2 and u 2 2 respectively. Their trajectories on 
the (x,t)—plane intersect near the point (X,T), which can be found from the system of equations 

(c.f. m) 


u\X + iq l T = — In 
1 2 


1 


I pi - q 1 

Pi + 9i 


V2X + u 2 1 T = - In 


| P2 ~ 92 

P2 + 92 


Far from this point as t —>• ±00 the solution tends to the sum of two simple kinks moving along the 
trajectories 

Vl (x — xf) + v^ l t = 0, V 2 ^X — x 2 ) + V 2 l t = 0 

respectively. Thus, the effect of the interaction is a shift of a straight line trajectory of each kink. 
In Figure [2 we plot a two-kink solutions when n = 2 for aq and a 2 at time t = 6 with 

v i = 0.9, (ai = (i,0,0,1); v 2 = 1.1, (a 2 = (i, 0,1,0). 

These two kinks intersect at the point x = 0. At t = 6, the position of the kink corresponding to 
v\ is at x =— —7.41 and the other kink is at x = —4.96. 

In order to determine the shift X 2 = x 2 — x 2 of the trajectory of the second kink we consider 
the limits t —>• ±00 assuming p 2 = V 2 X + v 2 1 t to be finite. In this case 


— 1 ^1 L I L/ \ \ 

pi = v\x + t = — p2 -\ -1-9 —>• ±00, as M ±00. 

Z^2 v\ ' 


t 


Therefore 


Ti(A, x,t)a 2 ) 


<h(A,iiq, T 0 (mi,x,t)ai))T 0 (A,x,t)a2) ->■ (A)T 0 (A, x, t)a 2 , as t —> ±00, ( 57 ) 


where 


$±(A) 


lim $(A, ivi, 'S>o(ivi,x, t)a±}). 

pi —>-±oo 
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The latter limit can be easily computed since t, t ) is a function of p\ only (1381) . In particular 

we have 

lim Q'Z>o(ivi,x,t)ai)(ai'S>o(ivi,x,t) T _ 1 
pi->±oo (aiT 0 (i^i, t, t) T Q^o(ii/i, x,t)ai) 2 

and therefore 


(l =R 0 \ 

I ±i 1 0 
\ 0 0 0 / 


$f(A) 


1 

A^TX 2 


A 2 — vf ±2Aiq 0 

=F2A^i A 2 - ul 0 

0 0 (A 2 + vl)I 


(58) 


Noting that (A), \ko(A,T,i)] = 0 we can represent the limits (1571) in the form 


^f(iu 2 )^ 0 (iu 2 ,x,t)a 2 ) = ^0 {iV 2 ,X,t)d 2 


±\ 


where 


a 0 ) = 


v\ + v\ 

7l 2 2 ’ 

U 2 ~ 


Using formula 


i{niP 2 T 72 92) 

7 i«2 =F I2P2 

&2 

we obtain 

x$ = lim xo = 

" p —»+00 

t 7 = lim To = 

p —> — OO 

Therefore, the position shift for the kink with speed 1 jv\ is 


2V\V2 

T2 2 2 

7/^ _ 7/^ 


1 In! 

(71 +72XP2 - 92), 

2 v 2 1 

1 In! 

(71 -72XP2 + 92) 1 

(71 72 ) (P 2 92 ) | 

2 v 2 1 

(71 + 72 ) (P 2 + 92 ) 


AT2 = xt ~ To = 


1 , / 71 + 72 \ 2 1 , / iq + ^2 \ 2 2 , v 2 + V\ 


■ In 


= — In 

2v 2 V7l —72/ ^2 V^l - y 2) V2 V 2 —V\ 

In a similar way one can obtain that the position shift for the kink with speed 1 jv\ is 


= — In ■ 


Ati = x± — x l = — — In f 1/1 U2 
V\ \ V2 ~ V\ 


2 v 2 -v\ 

= — In ■ 


vi v 2 + 

Notice that these position shifts Ati and At 2 for kinks after collision only depend on the pole 
position zq and v 2 . They are exactly the same as that of the scalar sine-Gordon equation ([3]) given 
in 


3.2 Interaction of two breathers 

A two-breather solution corresponding to the sets 

{Pi,ai) = {pi,qi,a\,... ,ai) T }, {p 2 ,a 2 } = (p 2 , q 2 , a \,..., a 2 ) T }, 0 < \m\ < \p 2 \ 

represents two breathers moving to the left with speeds |/ii| -2 and I/12I 2 respectively. Their 
trajectories on the (t, t)-plane intersect near the point (X, T), which can be found from the system 
of equations (c.f. (|55l) l 

uiX+ui^f + = \ \n\ Pl Qll \, 1 = 1,2. 

2 Pi+qii 
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Far from this point as t —>• ±oo the solution tends to the sum of two simple breathers moving along 
the trajectories 

* - xf + ill + )~H = 0, 1 = 1,2. 

with phase (jrf 1 . Thus, the effect of the interaction is a shift of a straight line trajectory and the 
change of the phase for each breather. 

In order to determine the shift 772 = V 2 ~ *?2~ °f th e trajectory and the change of the phase for 
the second breather we consider the limits t —» ±00 assuming P 2 = x + (y| + to be finite. 

In this case 


pi = x + (71 + w 2 ) = p 2 + 


7i + w? 


1 

72 + w 2 


t —>• ±00, 


as f —>• ±00. 


Therefore 


^i(A, x,t)a 2 ) = $(A,/ri,^o(/ri,x,t)oi))^ , o(A,x,t)a 2 ) -»• (A)^o(A,x,i)a 2 , as i —> ±00, (59) 

where 

(A) = lim $(A,//i,^ 0 (//i,x,i)ai)). 

pi ^-±00 

The latter limit can be easily computed using \&o(mi ,x,t) (15(71) . In particular we have 

lim A = ^ 

p i->-±oo ‘-fl 



and therefore 


<f>±(A) = 


/ |A 2 -|mi| 2 ) 2 — WUj , 4 Ao;i(A 2 -|mi| 2 ) n \ 

(A 2 —Mi) (A 2 —Mi) vw '° " v ' U 


V 


. _ . . .. '(A 2 -Ml) (A 2 -Ml) 

4A^i(A 2 -|mi| 2 ) (A 2 —Imu) 2 —4A 2 £a> 2 

+ U2-„2ia2 _,i-M (A 2 -m?)(A 2 -m 2 ) 

0 


(A 2 —M 1 KA 2 —Mf) 

0 


0 


(60) 


Noting that [$1 (A), 'h 0 (A,x,t)] = 0 we can represent the limits (15111) in the form 

$f(^ 2 )^o(V 2 ,x,t)a 2 ) = ^ 0(^2 ,x,t)af). 

We denote (1,1) and (1,2) entries of matrix 3 >i~(m2 ) by K\ and K2, respectively. Notice that 

(m 2 “ lMi| 2 “ 2/i 2 wi*) 2 {p2 ~ Mi) 2 (M2 + Pi) 2 


Ki — K2* = 


Kl + K 2 * = 


(m! — miXm! — m 2 ) (m!-m?)(m1-m?) ’ 

(M2 - lMl| 2 + 2/72W1*) 2 _ (M2 + Ml) 2 (M2 - Ml) 2 


(mI-m?)(mI-m?) 

Using formula for 770 in (1551) we obtain 
, 1 




Xn = 


lim 7/0 = 


iO 2 pi - H-OO 

Xn = — lim t/o = 

U >2 Ml —>—oo 


1 In 

(«i + K 2 i)(p2 - 92*) | 

- 1 In 

(mi + M2) 2 (mi - M 2 ) 2 ( 9*2 - 92*) | 

2w 2 1 

1 In 

Kl - K 2 i)(p 2 + 92*) 1 
(«1 - K 2 *)(p2 - 92*) | 

2w 2 1 

- 1 In 

(mi - M 2 ) 2 (mi + M 2 ) 2 (jP 2 + 92*) 
(mi - M 2 ) 2 (mi + A* 2 ) 2 (M 2 - 92*) | 

2w 2 

Kl + K 2 i)(p2 + 92*) ' 

2w 2 

(mi + M 2 ) 2 (mi - M 2 ) 2 (jP 2 + 92*) 
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Therefore, the position shift for the breather with speed 1 /l/J-21 2 is 


Ax2 = X2 — x 2 



(Mi + m)(m -m) 

(m -m)(£i +m) 


Using the formula for 00 in (1561) . we compute the phase shift for this breather. Indeed, we have 


<Pt 


0 2 


1 . , _ 1 ~ ^2) 2 (^l + m) 2 (P2 + 92^ . 

P 1-+00 ° ^ \(/Xi + M2) 2 (Al - H2) 2 (P2 ~ Q2i)J ’ 

1. , „ /" (m + ^2) 2 (Al ~ m) 2 (P2 + g2») \ 

P 1-+00 0 arg \(m -/r 2 ) 2 (m +/r 2 )2(p 2 _<? 2 *)J 


Hence the phase shift is 


A02 = 02 


02 


\(m + m) 2 (m -m) 2 / 


4 arg 


0 (m -m)(m + m) \ 
V(mi + m)(m - m)/ 


In a similar way we obtain that the position and phase shifts for the breather with speed l/|/xi| 2 
are 


Axi = x\ — x 1 
A0i = 0+ - 0“ 


2 lnl ^ 1 ~^(f il + f i 2 h . 
wi (m + m)(m - m) ’ 

1nr:r Am + m)(m-mA 

s \(m - m)(m + m)/ 


Similar to the case of interaction of two kinks, the position and phase shifts for breathers after 
collision only depend on the pole position and p, 2 , which are exactly the same as that of the 
scalar sine-Gordon equation Q given in M- 


3.3 Interaction of one kink and one breather 


A kink-breather or breather-kink solution corresponding to the sets 

{iv,a) = {ip,q,a},...,a n ) T }, {p, a) = (p, q, a 1 ,..., a n ) T }, 0 <v<\p\, 


where v,p, q, a 1 ,..., a n E R and // = 7 +iu,p, q, a 1 ,..., a n E C, represents one kink and one breather 
moving to the left with speeds and |/u| -2 respectively. Their trajectories on the (x,f)-plane 
intersect near the point (X,T), which can be found from the system of equations (cf. ([391) 1 


vX + v~ x T = 



P-q 

p + q 


uX + w(7 2 + w 2 


- 1 T = - In 
2 


p — qi 
p + qi 


Far from this point as t —>• ±00 the solution tends to the sum of one kink and one breather moving 
along the trajectories 


v(x — x ± ) + v l t = 0, (x — x^ + (7 2 + cj 2 ) 1 t = 0 

respectively, and the phase for the breather is 0 :t . Thus, the effect of the interaction is a shift of a 
straight line trajectory of each of them. 
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In order to determine the shift Ax = x + — x~ of the trajectory and phase shift A<0 = 0 + — 0“ 
of the breather we consider the limits t —y ±oo assuming p = x + (y 2 + u 2 )~ 1 t to be finite. In this 
case 

1 1 


p = vx + v = z; I p + I — — 


+ uj 1 


t —> ±oo, as t ±oo. 


Therefore 

'I'i(A, x, t)a) = <h(A, iu, x, t)a))To(A, x, t)a) —>• x, t)a, as t —>• ±oo, (61) 

where 


1 


A 2 — z / 2 ± 2 Az 9 

v 2 ,,2 


(A) = lim <f>(A, iv, 'I'oUz', x, t)a)) = —5 - 77 =p 2 A^ A 2 — v 

p-> ±00 A z + z/ z 1 


0 


0 (A 2 + v 2 )I 


according to (1581) . Noting that [^(A), To(A, x, t)] = 0 we represent the limits (fiTTIl in the form 


( 11 )^ 0 ( 11 , x,t)a) = T 0 (^,x,t)a ± ) 


where 


(a ± = 


p 2 - v 2 


p ± 


2pu 


9 1 9 /' 1 9. 9" 1 / 5 T" 9 . 9 f 9. 9 

;U Z + l/- /U z + Z/ Z /U + V A p A + V A 


2pu 


tP + 


P 2 -v 2 ~ -1 


9,a 


Using formula (|55D we obtain 


£+ = I lim 

uj p-H -00 2oj ((z/— a ;) 2 + 7 2 )|p + 

r~ = - lim n - 1 1n ((^ ~ ^) 2 + 7 2 )|p ~ g»| 

w p-y -00 ° 2 w ((z/+ a ;) 2 + 7 2 )|p + ’ 

It follows that the position shift for the breather is 


Ax = x+-iT = -ln |M + ^ 


1 {v + lo) 2 + 7 2 

= — In ■ 


w ' p — vV u ; (u — oj ) 2 + 7 2 


Using formula (|56l) we obtain 


= lim 0 o = arg 

p— H-oo 


(p — ui) 2 (p + qi) 

(/i + z/*) 2 (p - qi) 

Hence the phase shift for the breather is 

A0 = 0 + — 0 _ = 4 arg 


= lim 0 o = arg 

p—> — OO 


H~ vi 


(p + Z9*) 2 (p + qi) 
(p — vi) 2 {p — qi) J 


p + ui 


which are only dependent on the positions of poles, the values of v and p. 

In order to determine the shift Ax = x + — x~ of the trajectory of the kink we consider the 
limits t —» ±00 assuming p = vx + u~ 1 t to be hnite. In this case 


p = x + (y 2 + a; 2 ) = - - ( - 


1 


1 


i / 2 y 2 4- u 2 


=poo, as t-l ± 00 . 
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Therefore 


Ti(A, x, t)a) = T(A, //, To (fi, x, i)a))To(A, x, t)a) T±(A)To(A, x, t)a, as M Too, 
where 


/ 


(\ 2 -\ii\ 2 ) 2 -A\ 2 uj 2 


$1 (A) = lirn T(A ,n,^ 0 (n,x,t)a)) = 

p—>•=Foo 


<x 2 -ij: 2 )(x 2 -i i 2 ) T (\ 2 -^ 2 )(\ 2 ~p, 2 ) 

A \/.i7 \2_I,, |2\ 7 _I / ( 12\2_ a \ 2, .2 


4Aa;(A 2 —|/i| 2 ) 


± 


4Acj(A 2 —|p.| 2 ) 


V 


(X 2 -^i 2 ){X 2 -p, 2 ) 

0 


(A 2 —|/t| 2 ) 2 —4A 2 o; 2 
(. X 2 -p 2 )(X' 2 -fi 2 ) 

0 


(62) 


according to (1601) . Noting that [T^A), Tq(A, x, t)] = 0 we represent the limits (|62l) in the form 


^f(iix)^/ 0 (iu, x,t)a) = Tq (iis,x,t)a 


±\ 


where 


a ± ) = 


/ • ((7 2 +uJ 2 -t-i' 2 ) 2 +4i' 2 Lj 2 )p±4isLj('y 2 -\-u 2 +v 2 )q \ 

* ('y' 2 +oj 2 +v 2 ) 2 —4v 2 u 2 

±4i'o;(7 2 +o; 2 +^ 2 )p+ ( (•y 2 +u> 2 +u 2 ) 2 +4u 2 oj 2 ) q 
(i 2 +uj 2 +u 2 ) 2 -4u 2 u: 2 


Using formula (13911 we obtain 


-+ v i , ((y — w ) 2 + i 2 ) 2 \p - q\ 

p-+-o o 2v ((n + cu) 2 + 7 2 ) 2 |p + g| 

V 1 ! ((^ + w) 2 + 7 2 ) 2 |p — g| 

x = iim xn = — In —- -77 -tttttt -r • 

p->-l-oo 2 v {[v — to) 2 + 7 2 ) 2 \p + q\ 


Therefore, the position shift for the kink with speed 1 jv 2 is 


2\4 


Ax = x + — x = 


In 


((*-< + Y) 


= — In 


(v - to) 2 + 7 2 


2v ((n + w) 2 + 7 2 ) 4 n (u + cn) 2 + 7 2 


In Figure [3l we give the contour plot of a\ to in (x,t)-plane to demonstrate the phase and position 
shifts of one kink and one breather interaction, where the kink and breather are corresponding to 


J 2 = 1, Kink = (10^,7,4,5.916), n = 2 + 0.6*, (a bre ather = (*, 2,3i, 2.45). 

They intersect approximately at the point (—0.930,0,0626). After the collision, the position of 
kink gets shifted forward by 0.455 and the position of breather gets shifted backward by 0.759. 


4 Conclusion 

In this paper we have found and studied kink, breather and multi-soliton solutions of the vector 
sine-Gordon equation (jTJ). The soliton solutions obtained in this paper are of rank one and are 
generic. Actually a kink solution cannot be of rank more than one. For example, in order to 
construct a kink solution of rank two, one needs to have two vectors 

(oi = (ipi,gi,af), (a 2 = (ip2,q2,a%), Pi,P 2 ,qi,Q 2 € M, ai,a 2 £i" 
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satisfying the conditions (aiai) = 0, ( 0202 ) = 0, ( 0102 ) = 0. This implies that vectors (a\ and (02 
are linearly dependent. However, in the case of breathers the corresponding vectors are complex 
and therefore there is no such obstacle. Thus higher rank breather solutions may exist. 

Our results show that a fc-kink solution of the n-component (a E M n ) vector sine-Gordon 
equation (k < n ) is also a solution of ^-component sine-Gordon equation. Similarly, a general 
/c-breather solution of n-component equation JT]) {k < n/2) is a solution of 2 k component vector 
sine-Gordon equation. In particular, a general breather solution of the vector sine-Gordon equation 
cannot be obtained from any solution of the scalar sine-Gordon equation, but is a solution of a 
two-component vector sine-Gordon equation. 

The result of interactions (positional shifts and phase shifts) due to collisions in the mi-kink 
and m 2 -breather solution of the vector sine-Gordon equation m is exactly the same as in the case 
of the scalar sine-Gordon equation Q, which has been studied in M- It only depends on the 
positions of poles in the dressing matrices. This is significantly different from the case of the vector 
generalisation of nonlinear Schrodinger equation, where the result of interactions depends on the 
initial polarisation of the colliding solitons m- This can be explained by the difference in the 
spectral and symmetry properties of their Lax operators. These properties have been essentially 
used in the construction of multi-soliton solutions. 


Acknowledgements 

The paper is supported by the Leverhulme Trust, AVM’s EPSRC grant EP/I038675/1 and JPW’s 
EPSRC grant EP/I038659/1. All authors gratefully acknowledge the financial support. JPW would 
like to thank S.C. Anco for providing some references. 


References 

[1] K. Pohlmeyer and K.H. Rehren. Reduction of the two-dimensional O(n) nonlinear u-model. 
Journal of Mathematical Physics , 20(12):2628-2632, 1979. 

[2] V. E. Zakharov and A. V. Mikhailov. Relativistically invariant two-dimensional models of field 
theory which are integrable by means of the inverse scattering problem method. Zh. Eksper. 
Teoret. Fiz., 74(6): 1953-1973, 1978. 

[3] H. Eichenherr and K. Pohlmeyer. Lax pairs for certain generalizations of the sine-Gordon 
equation. Physics Letters B, 89(l):76-78, 1979. 

[4] Ioannis Bakas, Q-Han Park, and Hyun-Jong Shin. Lagrangian formulation of symmetric space 
sine-Gordon models. Physics Letters B, 372(12):45-52, 1996. 

[5] Jing Ping Wang. Generalized Hasimoto transformation and vector Sine-Gordon equation. In 
S. Abenda, G. Gaeta, and S. Walcher, editors, Symmetry and Perturbation Theory, SPT 2002. 
World Scientific, 2003. 

[6] J. A. Sanders and Jing Ping Wang. Integrable systems in n-dimensional Riemannian geometry. 
Moscow Mathematical Journal, 3(4): 1369-1393, 2003. 


19 


[7] S.C. Anco. Hamiltonian flows of curves in symmetric spaces G/SO(N ) and vector soliton 
equations of rnKdV and Sine-Gordon type. Symmetry, Integrability and Geometry: Methods 
and Applications, 2:044, 2006. 

[8] S. C. Anco and T. Wolf. Some symmetry classifications of hyperbolic vector evolution equa¬ 
tions. Journal of Nonlinear Mathematical Physics, 12(Supplement 1):13-31, 2005. 

[9] Vladimir G. Ivancevic and Tijana T. Ivancevic. Sine-Gordon solitons, kinks and breathers as 
physical models of nonlinear excitations in living cellular structures. Journal of Geometry and 
Symmetry in Physics, 31:1-56, 2013. 

[10] S.V. Manakov. On the theory of two-dimensional stationary self-focusing of electromagnetic 
waves. Sov. Phys. JETP, 38(2):248-253, 1974. 

[11] Allan P. Fordy and Peter P. Kulish. Nonlinear Schrodinger equations and simple Lie algebras. 
Communications in Mathematical Physics, 89(3):427-443, 1983. 

[12] V.E. Zakharov and A.B. Shabat. Integration of nonlinear equations of mathematical physics by 
the method of inverse scattering. II. Functional Analysis and Its Applications, 13(3): 166-174, 
1979. 

[13] A.V. Mikhailov. The reduction problem and the inverse scattering method. Phys. D, 3(1& 
2):73-117, 1981. 

[14] L.D. Faddeev and L.A. Takhajan. Hamiltonian Methods in the Theory of Solitons. Springer 
Verlag, Berlin, 1987. 

[15] Stephen C. Anco, Nestor Tchegoum Ngatat, and Mark Willoughby. Interaction properties 
of complex modified Korteweg-de Vries (mKdV) solitons. Physica D: Nonlinear Phenomena, 
240(17):1378-1394, 2011. 

[16] A.V. Mikhailov. Integrability of a two-dimensional generalization of the Toda chain. JETP 
Lett., 30(7):414-418, 1979. 

[17] A.V. Mikhailov. Reduction in integrable systems. The reduction group. JETP Lett., 32(2):187- 
192, 1980. 



Figure 1: One breather solution. The curve ( x, oq, Q 2 ) with /i = 2 + 0.25i is plotted at time t = 3. 
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Figure 2: Two kink solution. The curve with v\ — .9, V 2 — 1.1 is plotted at time t = 6. 
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Figure 3: The contour plot of a\ in plane: Interaction of a kink with a breather 
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